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ELECTRON  BEAM  PROPAGATION  THROUGH  A  MAGNETIC 
WIGGLER  WITH  RANDOM  FIELD  ERRORS 


I.  Introduction 

Electron  beam  propagation  through  magnetic  wigglers  has  become  a  topic  of  recent 
concern  primarily  due  to  its  relevance  to  free  electron  lasers  (FELs).  Interest  in  FELs 
has  become  widespread1  since  the  FEL  offers  both  a  tunable  source  of  radiation  from  the 
microwave  to  the  subvisible  range  as  well  as  the  capability  for  producing  intense  power 
levels.  The  quest  for  FELs  to  serve  as  high  power  single  pass  amplifiers  has  led  to  the 
design  of  long  wigglers  extending  for  hundreds  of  wiggler  wavelengths.  The  design  and 
analysis  of  such  FELs  typically  assumes  the  wiggler  field  to  be  adequately  described  by  its 
ideal  sinusoidal  form.  However,  the  intrinsic  impreeisions  winch  occur  in  the  fabrication 
and  assembly  of  wiggler  magnets  result  in  a  magnetic  field  which  deviates  from  the  ideal 
sinusoidal  form  by  some  small  error,  SB.  Typical  state-of-the-art  wiggler  construction2 
results  in  intrinsic  field  errors  of  SB/BW  >  0.1%,  where  Bw  is  the  peak  ideal  wiggler  mag¬ 
netic  field  on  axis.  These  intrinsic  field  errors  lead  to  detrimental  effects  which  generally 
increase  with  the  axial  length  of  the  wiggler.3-6  Hence,  for  long  wigglers  the  detrimental 
effects  of  field  errors  become  exceedingly  important.  These  effects,  if  left  uncorrected,  may 
destroy  the  FEL  interaction  and  lead  to  a  loss  in  radiation  gain. 

Physically,  the  detrimental  effects  of  wiggler  errors  may  be  understood  as  follows.  As 
the  electron  beam  propagates  through  the  wiggler  in  the  axial  z-direction,  the  electrons 
encounter  a  series  of  errors  in  the  transverse  magnetic  field  <5B_l,  which  are  assumed  to  be 
random.  The  beam  electrons  then  experience  a  series  of  random  transverse  (vz  x  £Bj_ ) 
forces,  where  v-  cz  c  is  the  axial  velocity  of  a  relativistic  electron.  This  series  of  random 
forces  results  in  a  random  walk  of  the  electron  beam  centroid,  causing  the  centroid  to 
deviate  from  the  wiggler  axis  by  some  amount  Sz(z).  Statistically  speaking  (as  will  be 
discussed  below),  the  rms  magnitude  of  transverse  beam  displacement  6z(r)  generally 
increases  as  a  function  of  the  axial  propagation  distance  z.  In  an  FEL,  the  magnitude 
of  the  beam  displacement  Sz(z)  may  become  large  enough  to  prohibit  optical  guiding'  of 
the  radiation  which  decouples  the  electron  beam  from  the  radiation  beam  (the  electron 
beam  centroid  no  longer  overlaps  the  radiation  beam  centroid),  thus  destroying  the  FEL 
interaction;  or  worse  yet,  the  electron  beam  may  deviate  sufficiently  far  off  axis  so  as  to 
hit  the  wall  of  the  wiggler  magnets.  Ideally,  to  avoid  such  detrimental  effects,  it  may  be 
desirable  to  keep  the  magnitude  of  the  beam  w’alk  off  less  than  the  radius  rj,  of  the  electron 
beam,  (6x|  <  tv 

Not  only  do  the  wiggler  field  errors  cause  the  electron  beam  to  tvalk  off  axis,  the 
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errors  also  cause  the  parallel  energy  of  the  electron  beam  to  deviate  from  i t s  ideal  value 
(the  value  in  the  absence  of  field  errors).6  As  the  field  errors  induce  a  transverse  beam 
motion  (and  a  transverse  energy)  through  the  beam  walk  off,  the  errors  subsequently 
alter  the  parallel  energy  of  the  beam,  since  static  magnetic  fields  do  not  alter  the  total 
beam  energy.  Statistically,  the  variance  of  the  deviation  in  the  parallel  beam  energy  £q||(r) 
generally  increases  as  the  propagation  distance  r  increases.  Here,  7||  =  (1—  vr/c2)-3^2  is  the 
relativistic  factor  associated  with  the  axial  electron  motion.  In  an  FEL,  the  parallel  energy 
deviation  67||(r)  may  increase  in  magnitude  to  the  point  where  the  wave-particle  resonance 
is  no  longer  maintained,  thus  destroying  the  FEL  interaction  and  reducing  the  overall  FEL 
gain.  In  order  for  the  electrons  to  maintain  their  resonance  with  the  radiation  field,  it  is 
necessary  for  the  electron  parallel  energy  deviation  to  be  small  compared  to  the  intrinsic 
power  efficiency  rj  for  FELs  in  the  low  or  high  gain  regimes,8  |^7||/7||ol  <  *7,  where  7||0  is 
the  electron  parallel  energy  in  the  absence  of  field  errors.  In  the  trapped  particle  regime, 
the  parallel  energy  deviation  must  be  small  compared  to  the  ponderomotive  potential 
created  by  <he  FEL  radiation:9  )^7||/7||o|  <  |e4>p/(7Tnoe2)|,  where  4>p  is  the  ponderomotive 
potential  and  mo  is  the  electron  rest  mass. 

The  effects  of  random  field  errors  in  magnetic  undulators  were  first  analyzed  by 
Kincaid.3  In  his  work,  Kincaid  was  concerned  with  how  these  errors  affected  the  spon¬ 
taneous  radiation  spectrum  resulting  from  the  passage  of  an  electron  beam  through  the 
undulator.  Kincaid  studied  the  transverse  orbit  of  a  single  electron  in  the  1-D  limit,  ne¬ 
glecting  the  effects  of  transverse  weak  focusing.  Kincaid  also  assumed  a  specific  model 
for  the  random  errors,  in  which  the  axial  dependence  of  the  field  error  associated  with 
a  given  pole  pair  was  assumed  to  be  a  sinusoid  extending  over  half  an  undulator  period 
(see  the  discussion  in  Sec.  VI).  Elliott  and  McVey  also  analyzed  the  effects  of  field  errors 
in  undulators  and  wigglers.4  Again,  they  were  primarily  concerned  with  how  these  errors 
affected  the  spontaneous  radiation  spectrum  (for  undulators)  or  the  radiation  gain  (for 
wigglers).  Elliott  and  McVey  presented  theoretical  calculations  of  the  transverse  orbit  of  a 
single  electron  including  the  effects  of  transverse  focusing  based  on  a  model  which  assumed 
that  the  electrons  received  discrete  independent  velocity  kicks  from  field  errors  at  each  pole 
pair.  These  theoretical  calculations  were  then  supported  by  the  results  of  an  FEL  simu¬ 
lation  code  in  which  the  electron  dynamic  equations  were  averaged  over  a  wiggler  period. 
Shay  and  Scharlemann  used  a  similar  FEL  simulation  code  to  study  the  effects  of  field 
errors  on  FEL  performance.5  The  random  walk  of  the  beam  including  3-D  focusing  effects 
was  briefly  discussed  and  plots  showing  how  field  errors  reduce  FEL  output  power  were 
presented.  The  works  of  Kincaid,  Elliott  and  McVey  as  well  as  Shay  and  Scharlemann  all 
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discussed  how  the  detrimental  effects  of  field  errors  could  he  reduced  bv  periodic  external 
steering  of  the  electron  beam.  None  of  these  works,  however,  directly  calculated  the  effects 
of  random  fieid  errors  on  the  parallel  energy  of  the  electrons,  nor  did  they  perform  particle 
simulations  of  beam  transport  in  which  the  electrons  are  modeled  by  the  full  relativistic 
Lorentz  equations  (as  opposed  to  a  spatial  averaged  version). 

In  the  following  sections,  a  comprehensive  theoretical  and  numerical  analysis  is  pre¬ 
sented  determining  the  effects  of  random  wiggler  field  errors  on  relativistic  electron  beam 
propagation.  Specifically,  how  random  field  errors  affect  the  transverse  position  of  the  elec¬ 
tron  beam  centroid  as  well  as  the  parallel  electron  beam  energy  will  be  determined.  This 
study  concerns  the  effects  of  random,  homogeneous  errors  in  the  transverse  magnetic  field 
SBx(z)  which  statistically  have  zero  mean,  a  finite  variance  and  a  finite  autocorrelation 
distance.  It  is  assumed  that  the  field  errors  SB ±(~)  are  a  function  only  of  the  propagation 
distance  z,  which  is  a  valid  approximation  provided  that  the  transverse  variations  in  the 
field  errors  are  small  over  the  transverse  spatial  extent  of  the  electron  beam.  However, 
3-D  effects  are  retained  in  the  beam  dynamics  and  in  the  form  for  the  ideal  wiggler  field 
which  enables  the  effects  of  transverse  weak  focusing  to  be  studied.  In  this  analysis,  it 
is  assumed  that  there  exists  an  ensemble  of  wiggler  magnets  which  contain  statistically 
identical  magnetic  fields.  Expressions  are  derived  for  an  electron  beam  quantity  Q  for 
a  particular  wiggler  realization  (a  particular  member  of  the  ensemble)  as  well  as  for  the 
appropriate  statistical  averages  (denoted  by  angular  brackets)  over  the  members  of  the  en¬ 
semble,  such  as  the  mean  (Q)  and  variance  squared  ( Q 2)  —  (Q)2 .  The  theoretical  analysis 
presented  below  assumes  that  the  relativistic  electron  beam  dynamics  may  be  adequately 
described  by  the  dynamics  of  a  single  electron  located  at  the  beam  centroid.  This  assump¬ 
tion  is  then  supported  by  performing  full  scale  nonlinear  3-D  simula’  ions10  of  electron 
beam  propagation  including  finite  emittance  and  space  charge  effects. 

The  remainder  of  this  paper  is  organized  as  follows.  Section  II  of  this  paper  presents 
an  analytical  treatment  of  the  basic  properties  of  general  random  field  errors.  In  Sec. 
Ill,  an  analytic  theory  is  developed  for  electron  motion  in  a  helical  wiggler  with  general 
random  errors  in  the  1-D  limit.  This  theory  is  generalized  to  include  the  3-D  effects  of 
transverse  weak  focusing  in  Sec.  IV.  In  Sec.  V,  analytical  results  are  presented  for  electron 
propagation  through  planar  wigglers  for  the  case  of  flat  pole  faces  as  well  as  for  the  case  of 
parabolic  pole  faces.  Section  VI  presents  the  results  of  the  beam  propagation  simulation 
code.  This  paper  then  concludes  with  a  summary  and  discussion  of  the  results  in  Sec.  VII. 


II.  Random  Field  Error  Statistics 


In  an  actual  wiggler  magnet,  the  magnetic  field  will  deviate  from  the  ideal  theoretical 
sinusoidal  form,  Bu.,  by  some  small  amount  8B.  Hence,  the  total  magnetic  field  will  be 
denoted  by 

B(x,y,z)  =  B u,(x,y,z)  -1-  6Bx(z)ex  +  8Bv{z)ey.  (1) 

Throughout  the  following,  the  dependence  of  8B  on  the  transverse  coordinates  will  be 
neglected.  This  assumes  that  ( 6B(rmax,z )  —  8B(0,  z))~/6B(0,  z)2  <  1,  where  rmai  is  the 
maximum  transverse  displacement  of  the  electron  beam  from  the  axis. 

In  the  following  analysis,  it  is  assumed  that  there  exists  an  ensemble  of  wiggler  magnets 
in  which  the  associated  field  errors  £B  exhibit  known  statistical  properties.  Specifically, 
it  is  assumed  that  (£B(z))  =  0  and  that  the  correlation  functions  for  the  field  errors, 
{ 8B(z)8B(z  +  A:)),  are  known.  Furthermore,  it  is  assumed  that  the  random  field  errors 
are  homogeneous,  that  is,  the  correlation  functions  { 6B(z)8B(z  +  Az))  are  only  a  function 
of  A z.  These  correlation  functions  are  assumed  to  exhibit  the  following  properties: 


+ A*))  j  = 

<«,(--)«?,(*  + A*))  {*£*?>■ 

(8Bx(z)8Bv(z  +  Az))  =  0. 


for  Az  =  0 
for  |Az|  >  zex , 

for  Az  =  0  (2) 

for  |Az|  >  zcy, 


In  the  above  expressions,  ( SB 2)  and  ( SB 2)  are  assumed  to  be  constants  and  represent  the 
mean-squared  field  errors  of  the  x  and  y  components  of  the  magnetic  field.  Also,  zex  and 
zcy  are  the  autocorrelation  distances  for  the  errors  8Bx(z)  and  6Bv(z),  respectively,  which 
are  defined  by  the  expressions 


<£Az  {6Bx(z)8Bx(z  +  Az)) 
<£Az  { 6By(z)6By(z  +  Az)) 


(SBl), 

5c 


(3) 


Physically,  zc  is  the  distance  over  which  the  field  error  8B(z)  remains  coherent  (see  Fig. 
1).  Typically,  one  expects  zc  ~  Au,/ 2,  where  \u,  is  the  wavelength  of  the  wiggler  field. 
This  implies  that  the  scale  length  for  the  field  error  associated  with  a  given  magnet  pole 
is  approximately  equal  to  the  width  of  that  pole.  The  last  expression  in  Eq.  (2)  indicates 
that  the  x  and  y  components  of  the  field  errors  are  assumed  to  be  uncorrelated. 
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Another  quantity  of  interest  is  the  vector  potential  ^.4r,v(r)  associated  with  t lie  field 
error  SBXtV{ r).  Defining  the  normalized  vector  potentials  ba  =  c6A/(moC2)  and  au,  = 
cBu, / (ku,m0c2),  where  kw  —  2ir/\u>  is  the  wiggler  wavenumber,  then 


6az  = 


Bx 


i; 


dz'  bBv{z' )  and  6 


ay  =  ~~B7 1 


dz'SBAz'). 


(4) 


It  then  follows  that  the  correlation  of  the  normalized  vector  potential  errors  is  given  by11,12 

alkj 


(6ax{zx)6at{z2))  =  ^  J  '  dz1  ’  dz"  ( 6Bv(z')bBv(z ")) 


.2  1,2 
"'ll1 

B2 

^  tn 


/*2w  AOC 

/  dz'  dAz  {SByiz^SB^z1  +  Az)) 

J  —  oo 


1  +  0  ^ 


(5) 


_  on  D  -1.* 

—  zJz  —  2  guj  —  *• 


i92 

hi 


u)-ev> 


provided  |zm|  >  zcy,  where  zm  is  the  smaller  of  Zj  and  z2.  Throughout  the  following,  terms 
of  order  0(zc/z)  will  be  neglected.  Similarly,  for  the  y-component  of  the  normalized  vector 
potential  errors,  one  has 


(^av(z1)iay(z2))  =  2  Dy  zmi 


M> 

Bl 


k 2 


lP2Cti 


(6) 


provided  |zm|  >  zez.  The  above  expressions  are  used  below  to  calculate  the  statistical 
behavior  of  the  random  walk  and  parallel  energy  deviation  of  the  electrons. 


III.  Propagation  Without  Transverse  Focusing 

As  a  first  step  in  determining  how  field  errors  effect  electron  beam  propagation  through 
a  magnetic  w’iggler,  the  transverse  gradients  associated  with  the  wiggler  field  will  be  ne¬ 
glected.  Such  an  approximation  is  valid  provided  the  displacement  of  the  electron  beam  off 
axis  is  much  smaller  than  a  wiggler  wavelength.  This  amounts  to  neglecting  the  transverse 
focusing  force  the  electron  beam  would  normal}}’  feel  as  it  moves  off  axis  into  a  region  of 
stronger  wiggler  strength.  In  the  absence  of  field  errors,  this  transverse  gradient  in  the 
wiggler  field  gives  rise  to  betatron  oscillations.8,13  As  is  shown  in  Sec.  IV,  neglecting  the 
transverse  focusing  forces  is  valid  provided  (2 zkp)2  <$<  1,  where  =  au,ku,/(\/ 2j)  is  the 
betatron  wavenumber. 

The  ideal  wiggler  field  is  assumed  to  be  helical,  and  in  the  1-D  limit,  is  given  by 

B„(i)  =  Bw  (cos  kwz  e*  4-  sin  kuiz  ev) , 

Bw  (7) 

A„(z)  =  — (c°s  kvz  et  -f  sin  kwz  e„) , 

Ay) 

where  Aw  is  the  vector  potential  associated  with  B„.  Generalization  of  the  results  below 
to  the  case  of  a  planar  wiggler  is  straightforward  and  is  discussed  in  Sec.  V. 

Neglecting  the  transverse  gradients,  the  problem  becomes  1-D,  and  the  motion  of  an 
electron  through  the  magnetic  wiggler  is  completely  described  by  conservation  of  perpen¬ 
dicular  canonical  momentum  and  conservation  of  energy, 

e 

Pj.  =  “Aj.  and  7  =  7||7x,  (8) 

where  pj_  =  7 mvx  and  7x  =  (1  +  p\  i™2  c2  •  Note  that  7  is  a  constant  of  the  motion 
since  there  is  no  applied  electric  field.  Writing  vj_  =  vu,  +  6\±,  where  vw  =  ca.w/i  is 

the  electron  wiggle  motion  in  the  ideal  wiggler  field  Bu,,  then  the  deviation  6vj_  from  this 

ideal  motion  is  given  by  6vj_  =  (c/7)6a_L-  Notice  that  (6vj_)  =  0  and  that  the  mean-square 
perpendicular  velocity  deviation  is  given  by 

c2 

(^’l)  =  Dxz,  for  |r|  >  zCZtV,  (9) 

where  D j_  =  Dx  +  Dv.  Hence,  the  mean-square  perpendicular  velocity  of  the  electrons 
increases  linearly  with  the  distance  traveled  through  the  wiggler  due  to  the  presence  of 
field  errors. 
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A  more  important  quantity,  however,  is  the  transverse  displacement  of  the  electron 
orbit  off  axis.  This  is  obtained  from  the  perpendicular  electron  velocity  by  the  relation 
v_l  =  v:dx±/dz.  Again,  writing  x_l  =  x„,  4  ^x^,  where  x„.  is  the  electron  transverse 
wiggle  orbit  in  the  ideal  wiggler  field  Bu,,  then  the  transverse  orbit  deviation  is  given  by 

<5xj.  =  -  /  dz'f> aj.,  (10) 

7  Jo 

where  the  approximation  vz  =  c  has  been  made.  Again,  (fox)  =  0,  and  the  mean-square 
transverse  displacement  is  given  by 

(ixx)  =  for  |r|  >  zCZtV.  (11) 

7 

The  above  equation  indicates  that  the  rms  value  of  the  transverse  displacement  due  to 
field  errors  increases  as  z1^2 .  This  expression  agrees  with  the  results  of  Kincaid3  in  which 
the  random  walk  in  1-D  was  calculated  for  a  specific  field  error  model  (see  Sec.  VI). 


The  results  for  the  mean-square  values  of  the  electron’s  perpendicular  velocity  and 
orbit  given  in  Eqs.  (9)  and  (11)  may  be  understood  through  the  following  physical  ar¬ 
guments.  The  x-component  of  the  equation  of  motion  for  the  electron  in  the  combined 
magnetic  field  is  given  by 

d 

'ym0—vz  =  eBxl>sinkwz  +  e8Bv.  (12) 

at 

The  first  term  on  the  right  of  the  above  equation  is  the  force  due  to  the  ideal  wiggler  field 
whereas  the  second  term  is  the  random  force  due  to  the  field  errors.  Hence,  the  error  8BV 
produces  random  velocity  kicks  8vz.  Such  a  process  is  described  statistically  as  “velocity- 
space  diffusion”  and  one  expects  the  mean-square  velocity  to  scale  as  (8v2)  ~  2 Dxz.  where 
Dz  is  the  diffusion  coefficient.  Similarly,  6x  =  J  dt8vz  and,  hence,  6x  is  a  “time-integrated 
diffusion  process”,  which  exhibits  the  generic  form  (8x2')  ~  (2/3 )Dxzi .  One  should  keep  in 
mind  that  the  above  relations  only  hold  for  long  “times”,  that  is,  for  {z/ zc)2  >>  1,  where 
zc  is  the  autocorrelation  distance  for  the  random  force  8BV. 

Another  quantity  which  is  of  interest  is  the  parallel  energy  of  the  electrons  7||.  The 
statistical  behavior  of  7||  is  calculated  from  the  relation  7||  =  where  7  is  a  constant 

of  the  motion  and  7x  =  (1  +  ')'20\)112.  Here  0  =  v/c  is  the  normalized  electron  velocity. 
Using  the  above  results,  one  can  write  0X  —  0X0  4  80 x-  where  0±o  is  the  normalized 
perpendicular  velocity  in  the  absence  of  field  errors.  Defining  Ixo  =  (1  4  720lo)l/:  and 
7||o  =  1/7x01  then  the  parallel  electron  energy  is  given  by 


7 II  =  7||o 


f  +  7jjo  {20xo  '  Ox  d"  &0 1) 


1 


(13) 


For  t lie  present  case  of  electron  propagation  in  the  1-D  limit,  fllv  —  au./7  and  bf31  = 
6«_l/7-  The  statistical  properties  of  the  parallel  energy  are  easily  calculated  bv  expanding 
the  expressions  for  7||  and  -> jj  to  include  all  terms  of  order  fin2x.  Specifically,  the  mean 
(7|l)  and  the  square  of  the  variance  (7jj)  -  (7|| )2  of  the  parallel  energy  are  given  by  the 


expressions 


<7||)  --IL°  =  — Z—  f  f 1  -  (Dz  +  Dv)  -  lai(Dt  -  Dy)cos2ku,z]  ,  (14a) 

7||o  7io  LV  2  /  2 


<7n)  -  <7||)' 


2 

=  (Dz  4-  Dy)  +  (Dt  -  Dv)cos2fcl 

~y 1 „ 


(146) 


where  7x0  =  (1  +  a^,)1^2.  In  the  limit  Dz  =  Dy  (( bB\)zcz  —  (6B2)zcy),  the  above 
expressions  reduce  to 


(7||)-7||o  (l  ~  fll/2) /<*£ 


(1  +  «= 


2  1  2  _ 

^  t<*  *“  ) 


(*r{)  -  <7|i)3 


(1  +«lf\Bi/  ■ 


(15  a) 


(156) 


As  before,  the  above  formulas  apply  provided  ( z/zc )2  >>  1.  Notice  that  (7(|)  —  7||0  may 
either  be  positive  or  negative,  depending  on  whether  or  not  a^/2  >  1.  This  indicates  that  it 
is  possible  for  the  field  errors  to  perturb  the  electrons  in  such  a  way  that  the  perpendicular 
energy  in  the  ideal  wiggle  motion  of  the  electrons  is  converted  into  parallel  energy.  Also, 
notice  that  the  variance  of  the  parallel  energy  deviation  increases  as  y/z. 


IV.  Propagation  with  Transverse  Focusing 


Transverse  focusing  of  an  electron  beam  results  from  the  transverse  gradients  in  the 
wiggler  magnetic  fields.  Specifically,  as  an  electron  beam  moves  off  axis  it  also  moves 
into  a  region  of  higher  magnetic  field  which  tends  to  focus  the  electron  beam  back  toward 
the  axis.  This  weak  focusing,  in  the  absence  of  field  errors,  produces  beam  oscillations 
at  the  betatron9’13  wavelength  =  \/2'y / aUp •  Hence,  in  the  presence  of  field  errors, 
one  expects  that  as  the  electron  beam  begins  to  random  walk  off  axis,  the  weak  focusing 
forces  tend  to  steer  the  beam  back  towards  the  axis,  thus,  diminishing  the  magnitude  of 
the  transverse  random  walk  of  the  beam  centroid. 

To  study  the  transverse  motion  of  an  electron  beam  in  the  presence  of  weak  focusing, 
a  helical  wiggler  field  is  assumed  with  a  normalized  vector  potential  given  by  the  following 
model  equations: 

ax  =  au  .(l  -f  k^y2  /2)  cos  kv,z  +  6ax(z ), 

(16) 

av  =  a„  (1  4-  k^x  / 2)  cos kwz  +  6ay(z), 

where  it  is  assumed  k2,x2  <<  1  and  k^y2  <<  1.  This  model  gives  an  adequate  approxi¬ 
mation  to  a  realizable  helical  wiggler  field  near  the  axis.14 

The  electron  motion  is  described  by  the  relativistic  Lorentz  equation,  which  may  be 
written  in  the  following  form: 


d  (  c  \  c  d 

Ti  v*  ~  jay  =  ~ 7 v  a 

d  f  c  \  c  d 

dt  V*  ~  7 a  /  ~~7V'5ya' 

The  above  transverse  equations  of  motion  will  be  solved  in  a  parameter  regime  assuming 
the  following  approximations.  Physically,  one  expects  the  perpendicular  electron  motion 
to  be  dominated  by  a  fast  wiggle  motion  t >±f  and  a  slower  random  motion  6v±  resulting 
from  the  finite  field  errors.  Hence,  it  will  be  assumed  6v  iA’i  /  <<  1.  Also,  cases  of 
physic^  ..  terest  occur  when  the  transverse  random  walk  6xj_  of  the  electrons  becomes 
much  r  er  than  the  amplitude  of  the  fast  wiggle  motion  Xj Hence,  it  will  be  assumed 
that  x j_,  <<  1.  As  will  be  shown  below,  these  two  inequalities  imply  that 


(17) 


<<  1. 
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(18) 


In  actuality,  the  results  derived  below  describing  the  random  motion  of  the  electrons  are 
somewhat  more  general  than  the  above  inequalities  imply.  However,  such  a  derivation 
becomes  too  detailed  to  be  presented  in  this  text.  Instead,  the  above  inequalities  will  he 
assumed,  which  corresponds  to  the  region  of  physical  interest,  which  greatly  simplifies  the 
derivation. 


The  x-component  of  the  electron  motion  is  calculated  bv  letting  vx  =  vTf  +  6vx,  where 
vZf  =  (ca„, /7)(1  -1-  k^,y7 /2)  cos  ku,z  is  the  fast  wiggle  oscillation  in  the  local  wiggler  field. 
Assuming  Sv\/v~{  <<  1  and  x^jSx1  «  1,  where  dSx/dt  —  6t>x  and  dxj/dt  =  vZf,  gives 

~;Sx  +  klSx  (1  —  cos  2 kwz)  =  —  -^-6ax,  (19) 

dz~  7  a: 

where  kp  =  awkw/{y/2 7)  is  the  betatron  wavenumber.  In  deriving  the  above  equation,  a 
change  of  variables  from  the  independent  variable  i  to  the  independent  variable  z  was  made 
along  with  d/dt  —  vzd/dz  ~  cd/dz.  Letting  Sx  =  Sx,  +  Sxf  and  assuming  Sx7/Sx]  <<  1, 
gives  Sx f  ~  —(kp/2kw)76x,  cos2kwz  and 


— Sx,  +  k76x,  = 


SJB, 
'  B„ 


(20  a) 


Clearly,  Sxj  represents  a  small  correction  to  Sx  on  the  fast  time  scale  d/dz  ~  2fcu,  which 
shall  be  neglected.  In  the  above  equation  for  Sx,,  the  second  term  on  the  left  represents 
the  weak  focusing  force  from  the  transverse  gradient  of  the  wiggler  field  and  the  term  on 
the  right  represents  the  random  force  from  the  finite  field  error. 


Similarly,  the  y-component  of  the  motion  can  be  calculated  by  letting  vy  =  vyf  +  />ry, 
where  vvj  —  {ca^J 7)(1  4-  k^x2  /2)  sinker,  and  assuming  Sv *  /•?/  <<  1  and  y'f/Sy 2  <<  1. 
Setting  Sy  =  Sy,  +  Syj  gives  Syf  =  (kp /2kw)2Sy,  cos  2 kwz  and 


d- 


6y,  +  k7fiSy ,  = - -k, 


SBX 


dz 2  "  '  7  Bw 

Solving  for  the  random  orbits  Sx,  and  Sy,  is  straightforward  and  one  finds 

Q'ixtk\p 


(20i>  l 


Sx,(z)  =  -■ 


Ifkp 


f  ~  j  1  •  1  ft  ^ 

/  dz  sin  kp  (2  —  z) — - - 

Jo 


ty 


J‘  dz’ smkg 


Bw 
fBx(z') 
B„  ' 


(21) 


Likewise,  the  random  perpendicular  velocities  Sj3x ,  —  bvx,jc  and  S0y,  =  Svy,/c  are  given 
by 

**,(=') 


6(3; 


f‘d={  costs  (=■---) 

1  Jo 


S0V.{Z)  =  - 


J  d~ 


COS  kp  ( r '  -  z) 


Bw  ’ 
fBz(:') 
Bu, 


(22) 


10 


Physically,  6t,  and  by,  represent  diffusing  betatron  orbits  which  are  the  result  of  random 
velocity  kicks  from  the  finite  field  errors  in  the  presence  of  weak  focusing  forces. 

In  order  to  calculate  the  mean-square  values  for  the  transverse  electron  orbits,  it  is 
necessary  to  evaluate  the  following  expression: 

/  dr'/  dz"  sin  kp  (r'  -  z)  sin  k0{z"  -  z)  (6B{z')6B{z" )) 

Jo  Jo 

~  /  dz'  f  dAz- [cos  kp  Ar  —  cos  kp  [2z'  —  2z  +  Az )] 

Jo  7—oc  2 

X  {6B(z')6B(z' +  Az))  [i  +  ^(t)]  (23a) 

-  J  dz'^[l  -  cos 2kp(z' -  z)]  (6 B7)  ze\l  +  O  +C(kpze)] 


where  (rc/r)2  <<  1  and  (kpzc)~  <<  1  has  been  assumed.  In  evaluating  the  above  expres¬ 
sions,  it  has  been  assumed  that  the  correlation  function  ( 6B{z)6B(z  +  Az))  is  independent 
of  z  and  becomes  zero  for  |Az|  >  zc,  as  is  discussed  in  Sec.  II.  Similarly,  one  can  show 


/  dz'  /  dz"coskp(z' -z)cosib/s(z"  ~z){6B{z')6B{z")) 

.  (236) 

These  expressions  are  then  used  to  calculate  the  statistical  averages  of  the  transverse  orbits 
and,  in  doing  so,  terms  of  order  0(ze/z)  and  0(kpzc)  will  be  neglected. 

Statistically  averaging  over  the  wiggler  ensemble  gives  the  following  expressions  for 
the  mean-square  quantities: 


<*>  -  o^y 


/*J£\h i 

\Bl/  2 


(sin2kgz\ 

Z~~2 kTJ' 


2  kp 


(24a) 

(246) 

(25a) 

(256) 
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Notice  that  the  above  expressions  reduce  to  the  corresponding  ]-D  expressions  in  the 

limit  (2kpz)2  <<  1.  The  above  expressions  indicate  that,  in  the  large  z  limit,  (^r*)^2  ~ 

~  z5^2  (as  opposed  to  z3^2  in  the  3-D  limit ).  This  argees  with  (lie  results  of  Elliott 

and  McVey4  and  of  Shay  and  Scharlemann.5  Hence,  weak  focusing  significantly  reduces 

the  magnitude  of  the  transverse  spatial  random  walk  of  the  beam  centroid.  Furthermore, 

notice  that  (6x2)  ~  {^Vl)  ~  1  /kp  and,  hence,  additional  external  focusing  (which 

increases  the  effective  value  of  kp)  subsequently  reduces  the  transverse  displacement  of  the 

beam  centroid.  Notice,  however,  that  {£/?*, )^2  ~  (S/37sY^‘  ~  z1/2,  as  is  the  case  in  the 

*>  \  f  *2. 

l-D  limit.  Hence,  in  the  large  z  limit,  weak  focusing  only  reduces  the  value  of  {Sfi7,) 

j  #2 

and  (f>fiys)  by  a  factor  of  root  two  as  compared  to  the  l-D  values. 

It  is  also  of  interest  to  determine  how  the  finite  field  errors  affect  the  parallel  energy 
of  the  electrons  in  the  presence  of  weak  focusing.  An  expression  for  the  parallel  energy  7|| 
is  easily  obtained  through  conservation  of  energy  7||  =  7/7 j_  provided  the  perpendicular 
motion  is  known,  7^  =  1  -f  72/?^.  Using  the  above  results,  one  has  fix  =  fix/  +  bfix»> 
where  fixj  is  the  fast  wiggle  oscillation  in  the  local  wiggler  field  and  6  fix,  is  the  random 
component  of  the  orbit  due  to  finite  field  errors.  It  is  then  straightforward  to  calculate  the 
various  statistical  moments  of  7||,  such  as  the  mean  (y\\)  and  the  square  of  the  variance 
(7jj)  —  (7||  )2.  This  is  done  by  expanding  the  corresponding  expressions  for  7||  and  7jj  to 
second  order  in  6B/BW.  One  finds 


('y\\)  — y»o 

7||o 


<zjj)  -  w 


3a2,  sin  2 k0zl 


sin  2  kpz 
2kp 


=  [{D*  +  D,)  +  (D,  -  I?y)cos2fcu.z]  (z  +  )  . 


cos2fct»z|,  (26a) 

(2  66) 


For  the  case  Dx  =  Dy,  the  above  equations  reduce  to 
(7||>  -7||o  7  -  ’  '  2 


7||o 


-(r&)  ©-[-(' *=?)•* 


3a2,  sin  2 kp-. 


2k  p 


(27a) 


<7||)  -  <7||)2  , 

(  Q-v'k\r 

I  '/SB'-Xaiz 

(  i  sin 2 kpz\ 

7no 

K 1  +  ) 

'  \Bl)  2  1 

l'+  2k,  ) 

(276) 


Notice  that  the  above  expressions  reduce  to  the  l-D  expressions  in  the  limit  (2kpz)2  <<  1. 
Furthermore,  in  the  large  z  limit,  both  (711)  ~  z  and  (7?)  —  (7||)2  ~  z,  as  was  the  case  in 
the  l-D  limit.  For  the  3-D  expressions  in  the  large  z  limit,  however,  (7||)  —  7||0  <  0  and 
(7jj)  —  (7||)2  has  been  reduced  by  a  factor  of  two  as  compared  to  the  l-D  value. 
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V.  Propagation  Through  Planar  Wigglers 

The  results  presented  above  are  for  wigglers  with  helical  fields.  It  is  straightforward 
to  generalize  these  results  to  planar  wigglers  with  linearly  polarized  fields.  The  electron 
motion  is  analyzed  using  the  methods  presented  in  the  previous  sections,  hence,  the  details 
of  such  calculations  will  not  be  repeated.  Below,  results  are  presented  for  planar  wigglers 
of  two  types:  i)  planar  wigglers  with  fiat  pole  faces  and  ii)  planar  wigglers  with  parabolic 
pole  faces. 

i)  Flat  pole  faces.  Consider  a  planar  wiggler  with  fiat  pole  faces  with  a  magnetic  field 
described  by  the  normalized  vector  potential 


a  =  aw  cosh  kwy  cos  kwzez  +  6az(z)ez  4-  Sav(z)ev. 


(28) 


The  above  vector  potential  gives  a  magnetic  field  primarily  in  the  y-direction  which  exhibits 
transverse  gradients  also  in  the  y-direction.  Hence,  intuition  indicates  that  the  electrons 
will  experience  wiggle  oscillations  in  the  x-direction  and  weak  focusing  in  the  y-direction. 

The  x-component  of  the  electron  motion  consists  of  fast  wiggle  oscillations  in  the  local 
wiggler  field  plus  random  velocity  kicks  in  the  absence  of  weak  focusing:  0X  =  0zj  +  60z, 
where  0zf  =  (aw/j)(l  +  k^y2 /2)  cos  k^z  and  60  =  6az/'y.  Hence,  the  random  part  of  the 
orbit  6x  is  described  by 


6x(z) 


and 


(29) 


The  y-component  of  the  electron  motion  consists  of  a  random  orbit  6y  including  the  effects 
of  weak  focusing.  One  finds  6y  =  6yt  +  6yj,  where  6yj  =  ( kp/4kl,)6y,  cos 2 kwz  is  a  small 
correction  to  6y  which  shall  be  neglected.  The  random  orbit  by,  is  described  by 


c  t  \  a.VfkyB  t  ,  i  ...  |  .  6Bz(z  ) 

6y,(z)  =  ——  /  dz  sin  kp  [z  -  z)  , 

lb  Jo  Bw 

/c  2\  _  ( avb  \  /  \  -cr  /  _  sin2kgr  \ 

'  W  \lb  )  \Bi  /  2  V  2b  }' 


It  is  also  straightforward  to  calculate  how  the  field  errors  affect  the  parallel  energy  of 
the  electrons.  One  finds,  to  second  order  in  \6B / Bw\,  that  the  mean  and  the  square  of  the 
variance  of  -y-jj  are  given  by 
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(31a) 

(316) 


where  7x0  =  (l+a«.  cos2  kvz)ll7 .  Notice  that  the  above  expressions  have  not  been  averaged 
over  a  wiggler  period. 


Hence,  for  flat  pole  faces,  in  the  x-direction  in  which  there  are  no  focusing  forces,  the 
rms  value  for  the  orbit  displacement  scales  as  z 3/2.  In  the  y-direction,  the  focusing  forces 
inhibit  the  orbit  walk  off  which  then  scales  as  z1^7  in  the  large  z  limit.  Notice  that  both 
the  mean  and  the  square  of  the  variance  of  the  parallel  energy  scale  as  z,  as  was  true  for 
the  case  of  a  helical  wiggler,  with  or  without  transverse  focusing. 


ii )  Parabolic  pole  faces.  Consider  a  planar  wiggler  with  parabolic  pole  faces  where  the 
normalized  vector  potential  is  given  by 


o-z 

av 


—  aw  cosh  (kwx/y/2?J  cosh  {k^y/y/^j  cos  kwz  +  6ax(z), 
=  ~aw  sinh  (kvx/y/2j  sinh  (kwy/ y/2^  cos  kwz  +  6ay(z). 


(32) 


Notice  that  this  is  for  the  special  case  of  equal  focusing  in  the  x-  and  y-directions;  that  is, 
kl  =  k7  where  k^  +  k7  =  k^,.  The  above  vector  potential  gives  a  magnetic  field  primarily  in 
the  y-direction  which  exhibits  transverse  gradients  in  both  the  x-  and  y-directions.  This 
indicates  that  the  electrons  will  experience  wiggle  oscillations  in  the  x-direction  and  the 
random  transverse  orbits  will  be  modified  to  include  weak  focusing  in  both  the  x-  and 
y-directions. 

The  x-cornponent  of  the  electron  motion  consists  of  fast  wiggle  oscillations  in  the  local 
wiggler  field  plus  random  velocity  kicks  in  the  presence  of  weak  focusing:  4-  f>l3x, 

where  ~  (a„,/7)(l  4-  fc*,r2/ 4  4  k7,y7 / 4)  cos  kvz.  The  random  conponent  of  the  orbit 
bx  can  be  written  as  6x  =  bx,  4-  bxj,  where  bxf  is  a  small  correction  to  the  random  orbit 
hxj  ~  (k^/Sk^bx,  cos  2 kvz,  and  where  bx,  is  described  by 
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Similarly,  the  y -component  of  the  motion  is  described  by  /?„  =  fiyj  +  6/3y,  where  (3yf  ~ 
— (o-u, /2y  )k^,xy  cos  k^z.  The  random  component  of  the  orbit  6y  can  be  written  as  6y  = 
8y,-\-8yjy  where  Syf  is  a  small  correction  to  the  random  orbit  6yj  ~  {k^/Sk^Sy,  cos  2 ku,z, 
and  where  6yt  is  described  by 


(34) 


Notice  that  the  above  expressions  for  the  random  orbits  6xa  and  Sy,  are  identical  to  those 
expressions  for  a  helical  wiggler  with  transverse  focusing  except  that  k@  now  must  be 
replaced  by  kp/y/2. 


It  is  also  straightforward  to  calculate  how  the  field  errors  effect  the  parallel  energy  of 
the  electrons  in  a  planar  wiggler  with  parabolic  pole  faces.  One  finds,  to  second  order  in 
\8Bf  J?u,|,  that  the  mean  and  the  square  of  the  variance  of  -yjj  are  given  by 


<7||)  ~  7110 
7|lo 


-  (v£,r.)  { (( if*' ) + (  bi  )=")  [-■ + 1 c°s 2Kz (2 

- -^)1  -  (kY-t  (*  -  i)  H w))- (35o) 


y/2kt 


<7(| )  ~  (7l| )2  =  /6B2y\Zl 

7j|o  V  7J.0  )  \Bl)  : 

where  7j.o  =  (1  +  a2,  cos2  fc^z)1/2. 


l1'^)  (2  + 


y/2kp 

sin  y/2k$z  ^ 
y/2kp  . 


(356) 


Hence,  for  parabolic  pole  faces,  transverse  weak  focusing  exists  in  both  the  x-  and 
y-directions  and  the  mis  value  of  the  orbit  displacement  in  both  transverse  directions 
scales  as  z1^2  (in  the  large  z  limit).  Notice  that  the  mean  and  square  of  the  variance  of 
the  parallel  energy  scale  as  z,  as  is  the  case  for  all  the  wigglers  examined. 
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VI.  Numerical  Simulations 

In  order  to  verify  the  above  analytical  theory  and  results,  a  numerical  code  was 
developed10  to  perform  full  scale  simulations  of  electron  beam  transport  through  mag¬ 
netic  wigglers.  This  code  is  a  fully  3-D  particle  simulation  which  includes  finite  beam 
eniittance  and  space  charge  effects.  This  code  simulates  beam  transport  in  either  helical 
or  planar  wiggler  configurations,  with  or  without  transverse  weak  focusing,  and  has  the 
capabilities  of  including  finite  wiggler  field  errors.  Once  the  magnetic  field  configuration  is 
specified,  electron  motion  is  simulated  by  solving  the  relativistic  Lorentz  force  equation. 

As  an  example,  electron  beam  transport  was  studied  in  a  helical  wiggler  with  finite 
field  errors.  The  magnetic  field  was  modeled  by  the  following  expressions: 


Bx  =  Bv  fid — ~~  J  cos kwz  +  -j1  (a:2  cos  kwz  +  xysh\kwz)  +  6Bx(z), 
\  (  k2  r2\  k2 

By  =  Bu,  l  1  + J  sin kwz +  (xy  cos  kwz +  y2  sin  kwz)  +  8By(z), 


Bz  =  Bwkw  ( 1  + 


(y  cos  kwz  —  x  sin  kwz) , 


(36a) 

(36fc) 

(36c) 


where  it  has  been  assumed  k^r2  =  fc^(  *2+y2)«l. 


The  functional  form  of  the  field  errors  8B(z)  was  chosen  as  follows.  It  is  assumed  that 
the  field  error  SBn(zn )  at  the  center  (x  =2„)  of  the  nth  pair  of  magnet  poles  is  a  random 
quantity  uncorrelated  with  the  field  error  8Bm(zm)  associated  with  the  center  of  the  mth 
pair  of  poles,  where  m  ^  n.  Furthermore,  in  order  to  preserve  the  continuity  of  8B(z)  as 
a  function  of  z,  it  is  assumed  that  the  field  error  6Bn(z)  associated  with  the  nth  pole  pair 
extends  over  the  region  zn  —  Xw/4  <  z  <  zn  +  A„,/4,  such  that  |62?n(x)|  is  maximum  at 
2  =  zn  and  is  zero  at  z  =  zn  ±  A.,/4.  For  simplicity,  the  following  functional  forms  axe 


Here,  A Bz  and  A Bv  are  the  maximum  field  errors  for  the  wiggler  in  the  x-  and  y- 
directions;  and  the  centers  of  the  pole  pairs  in  the  x-  and  y-directions,  znx  and  znv, 
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are  given  by  =  (n  -  l)A„./2  and  Cnj,  =  (n  -  l/2)Au./2.  Also,  c,lx  (or  e„y )  is  a  ran¬ 
dom  number  between  —1  and  1  which  is  constant  over  the  region  |r  —  rnt)  <  Au,/4  (or 
|z  —  :ny\  <  \v[4)  and  is  uncorrelated  with  the  value  of  emT  (or  Cn.y)  for  m  •/  «.  Further¬ 
more,  it  is  assumed  that  the  statistical  distributions  for  en*  (or  e„y)  are  identical  for  each 
pole  pair  (all  n).  The  distribution  for  £n  is  chosen  to  be  uniform  from  —1  to  1  such  that 
(e2t)  —  (fny)  =  1/3,  although  any  random  distribution  with  zero  mean  and  finite  variance 
would  be  equally  satisfactory.  This  model  for  the  field  errors  SBy  is  essentially  identical  to 
that  used  by  Kincaid,8  only  Kincaid  assumed  the  random  coefficients  eny  to  be  Gaussian 
distributed.  A  schematic  of  this  field  error  model  is  shown  in  Fig.  2. 

To  compare  the  above  numerical  model  for  6£(z)  to  the  analytical  model  which  as¬ 
sumes  6B(z)  to  be  a  random,  homogeneous  variable  with  autocorrelation  length  zc ,  it 
is  necessary  to  calculate  the  correlation  function  (6B(z)6B(z  -f  Az))  for  the  numerical 
model.  Recall  that  the  analytical  theory  assumed  SB(z)  to  be  homogeneous  such  that 
(SB(z)SB(z  +  Az))  is  a  function  only  of  A z.  For  the  numerical  model,  however,  this  is 
not  the  case  since  6B(z)  was  chosen  to  have  a  cos  kvz  (or  sin  kxz)  dependence  over  a  given 
pole  pair.  Hence,  a  comparison  of  the  numerical  model  with  the  analytic  theory  requires 
that  the  expression  for  the  numerical  model  of  ( 6B(z)£B(z  +  Ar))  be  spatially  averaged 
over  a  wiggler  period.  Doing  this,  one  finds  for  the  numerical  model,3 

(SB'(z)6B,(z  +  Az))  =  \aB\  <4)  (l  -  cos  kzAz 

11  A  <38) 

H — sinfcwlArl  ,  for  |Az|  <  ~, 

7T  2 

and  zero  otherwise  (along  with  the  corresponding  expression  for  the  field  error  correlation 
in  the  y-direction).  Using  the  definition  for  the  autocorrelation  length  ze  given  in  Sec.  II, 
one  identifies  that  for  the  numerical  model 


(iBl)  =  \ABl  (4>  =  jAB’, 

(IBl)  =  I A B\  (4)  =  iABj, 


(39a) 


(396) 


—  -  _  ___  \ 
-cr  —  ~cy  —  •)  Av. 

7T* 


(39c) 


In  the  above  expressions  for  ( SB 2),  the  factors  of  1/2  arise  from  the  spatial  average  of 
cos2  kxz  and  sin2  kwz.  Hence,  the  numerical  model  gives  an  rms  field  error  of  (SB2)1^2  = 
A B l y/6,  where  A B  is  the  maximum  field  error  of  the  wiggler  ensemble,  and  an  autocor¬ 
relation  length  equal  to  4\w/r2. 


The  above  model  for  the  wiggler  errors  was  used  in  t lie  numerical  code  to  simulate 
electron  beam  propagation  in  a  helical  wiggler.  A  particular  realization  of  a  wiggler  was 
obtained  for  a  given  single  set  of  random  field  error  parameters  cnx  and  e^y.  For  snrli 
an  individual  realization,  various  properties  of  the  electron  beam  were  calculated  such 
as  the  transverse  displacement  of  the  beam  centroid  and  the  variation  of  the  electron 
beam  parallel  energy.  These  beam  quantities  were  obtained  by  averaging  the  appropriate 
quantities  for  the  individual  electrons  over  the  distribution  of  electrons  within  the  beam. 
Ensemble  averages  of  a  beam  quantity  were  then  obtained  by  averaging  the  beam  quantity 
over  40  individual  wiggler  realizations  (40  sets  of  800  distinct  random  field  error  parameters 
tnx  and  e„y).  The  runs  described  below  are  for  a  helical  wiggler  of  length  L  =  40  m,  with 
Bw  =  4.3  kG,  Xw  =  5.0  cm,  and  (SB2)1^2  /Bw  =  0.3%;  and  for  an  electron  beam  of  energy 
7  =  270,  with  radius  rj,  =  0.08  cm  and  a  normalized  emittance  of  11.3  mrad-cm,  which 
matches  the  acceptance  of  the  transport  channel.  These  parameters  give  aw  =  2  and 
kp  =  9.5  m. 

A  typical  transverse  orbit  Sx(z)  occurring  in  a  single  wiggler  realization  is  shown  for  a 
case  without  transverse  focusing  in  Fig.  3  and  for  a  case  with  transverse  focusing  in  Fig.  4. 
Figure  3  shows  a  large  1-D  orbit  displacement  of  6.7  cm,  whereas  transverse  focusing  leads 
to  oscillations  about  the  axis  at  the  betatron  wavelength  with  a  maximum  displacement 
of  0.45  cm,  as  is  shown  in  Fig.  4.  Similarly,  the  parallel  energy  7||(z)  for  a  typical  wiggler 
realization  is  shown  for  a  case  without  transverse  focusing  in  Fig.  5  and  for  a  case  with 
transverse  focusing  in  Fig.  6.  Figure  5  shows,  in  the  1-D  limit,  a  2.5%  increase  in  the 
parallel  beam  energy.  This  is  in  agreement  with  Eq.  (15a)  which  indicates  an  increase 
in  the  mean  parallel  energy  provided  a2,  >  2.  Figure  6  shows  a  decrease  in  the  parallel 
energy  of  about  5%,  which  is  in  agreement  with  Eq.  (15b)  which  indicates  that  transverse 
focusing  will  lead  to  a  decrease  in  the  mean  parallel  energy.  Notice  also  that  the  parallel 
energy  in  Fig.  6  exhibits  oscillations  at  1/2  the  betatron  wavelength.  This  reflects  the  fact 
that  7||  depends  on  terms  proportional  to  6(3±  as  indicated  by  Eq.  (13). 

The  rms  transverse  displacement  of  the  beam  centroid  (6x2(z))1^2  as  obtained  from 
performing  the  ensemble  average  numerically  is  shown  in  Fig.  7  (solid  curve),  along  with 
the  theoretical  result  (dashed  curve),  for  the  case  without  transverse  focusing.  The  theory 
and  simulation  are  in  good  agreement  showing  a  maximum  rms  displacement  of  about  5 
cm.  Figure  8  plots  the  rms  beam  centroid  displacement  for  the  case  including  transverse  fo¬ 
cusing,  showing  the  comparison  between  the  numerically  obtained  ensemble  average  (solid 
curve)  and  the  theoretical  result  (dashed  curve).  Both  curves  exhibit  oscillations  at  1/2 
the  betatron  wavelength  and  show  a  maximum  rms  displacement  of  about  0.25  cm. 
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The  remaining  plots  describe  the  parallel  energy  ->||  of  1  lie  electron  beam.  The  nor¬ 
malized  ensemble  averaged  energy  (")'|| ) / Tpo  —  1  is  shown  in  Figs.  9  and  10  for  the  cases 
without  and  with  transverse  focusitig,  respecti vely,  comparing  the  .simulation  (solid  curve) 
and  the  analytical  theory  (dashed  curve).  Notice  that  in  the  absence  of  focusing  Fig.  9 
shows  an  increase  (of  about  1.3%)  in  the  mean  parallel  energy  whereas  Fig.  10,  which 
includes  focusing,  shows  a  decrease  (of  about  3.5%)  in  the  mean  parallel  energy.  This 
is  in  agreement  with  the  theory  for  the  case  au,  =  2.  Similarly,  the  normalized  variance 
<T|I  =  ((T(|)  ~  {'l'||)2)1^2/7||o  >s  shown  in  Figs.  11  and  12  for  the  cases  without  and  with 
transverse  focusing,  respectively,  comparing  the  simulation  (solid  curve)  and  the  analyti¬ 
cal  theory  (dashed  curve).  In  both  Figs.  11  and  12  the  variance  in  the  parallel  energy  is 
quite  large,  becoming  greater  than  10%  in  less  than  20  m. 

These  simulation  results  indicate  that  the  analytical  theory  gives  a  good  approxima¬ 
tion  to  the  qualitative  and  quantitative  behavior  of  the  electron  beam.  Several  differences 
which  appear  between  the  theory  and  the  simulations  may  be  attributed  to  the  following 
physical  and  l.americal  effects.  The  electron  beam  in  the  simulation  has  a  finite  emittance 
and  beam  cross-section.  This  results  in  a  mean-square  transverse  position  which  is  non¬ 
zero  initially.  This  is  seen  in  Figs.  3,  4,  7  and  8.  For  cases  in  which  weak  focusing  is 
included,  betatron  oscillations  occur  in  the  beam  motion.  Since  the  simulations  include 
finite  beam  emittance,  one  expects  the  amplitude  of  these  oscillations  to  be  larger  than 
the  predictions  ox  the  analytic  theory.  This  is  the  case  in  Figs.  8,  10  and  12.  Also,  the 
amplitude  of  this  oscillation  in  the  simulations  is  not  constant  as  the  beam  propagates  . 
This  effect  is  due  to  the  finite  sampling  size  used  in  performing  the  ensemble  average.  In 
principle,  an  infinite  set  of  wigglers  is  needed  in  order  to  replicate  the  smooth  functional  de¬ 
pendence  of  the  theoretical  ensemble  average.  For  the  ensemble  average  of  variables  whose 
variance  is  large,  the  sampling  size  used  must  be  large  in  order  to  recover  the  smooth  func¬ 
tional  dependence.  This  is  seen  in  the  simulations  shown  in  Figs.  9-12,  in  which  various 
statistical  moments  of  the  parallel  beam  energy  are  plotted.  These  ensemble  averages  were 
performed  over  a  set  of  40  realizations,  hence,  these  plots  do  not  exhibit  smooth  behavior. 
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VII.  Discussion 


The  above  analysis  indicates  that  intrinsic  magnetic  field  errors  perturb  electron  beam 
propagation  through  the  wiggler  and  lead  to  various  undesirable  effects  such  as  a  random 
walk  of  the  electron  beam  centroid3-6  as  well  as  fluctuations  in  the  parallel  electron  energy 
of  the  beam.®  These  detrimental  effects  were  studied  both  analytically  and  numerically. 
The  analytical  treatment  assumed  that  the  motion  of  the  electron  beam  could  be  ade¬ 
quately  approximated  by  a  single  particle  located  at  the  beam  centroid.  Expressions  were 
derived  for  the  centroid  motion  for  a  single  wiggler  realization  (a  particular  .occurrence  of 
random  errors  SB)  as  well  as  for  averages  taken  over  an  ensemble  of  wigglers  having  the 
same  statistical  properties.  The  field  errors  SB  were  assumed  to  be  random,  homogeneous 
functions  of  the  propagation  distance  z  with  zero  mean  and  with  a  finite  variance.  The 
transverse  dependence  of  SB  was  neglected  assuming  that  the  variation  in  SB  was  small 
over  the  transverse  extent  of  the  electron  beam.  Beam  propagation  was  studied  through 
both  helical  and  planar  wiggler  configurations,  with  and  without  transverse  focusing.  The 
results  of  the  analytical  theory  were  then  supported  by  performing  3-D  particle  simulations 
of  electron  beam  transport  including  the  effects  of  finite  emittance  and  space  charge.  The 
results  of  these  simulations  showed  good  agreement  with  the  analytic  theory. 

Physically,  as  the  electrons  propagated  through  the  wiggler,  they  experienced  ran¬ 
dom  \z  x  £B  forces  which  led  to  a  random  transverse  walk  of  the  beam  centroid.  In 
the  absence  of  transverse  focousing,  the  rms  value  of  the  centroid  displacement  scale  as 
{Sx2)1^2  ~  z3'2,  which  is  in  agreement  with  the  results  of  Kincaid3.  More  specifically, 
the  transverse  displacement  for  a  single  wiggler  realization  is  gWen  by  Eq.  (10)  and  the 
mean-square  value  averaged  over  the  wiggler  ensemble  is  given  by  Eq.  (11).  Transverse 
focusing,  however,  impedes  the  random  walk  such  that  ( Sx 2)*/2  ~  z1^2  in  the  large  z  limit, 
(2 kfiz)2  >>  1,  which  is  in  agreement  with  the  results  of  Elliott  and  McVey4  as  well  as 
Shay  and  Scharlemann.5  Expressions  for  a  single  realization  and  for  the  ensemble  averaged 
mean-square  of  the  centroid  displacement  with  transverse  focusing  are  given  by  Eqs.  (21) 
and  (24).  In  the  limit  (2 fc^z)2  >>  1,  the  rms  centroid  displacement  can  be  written  as 
{Sx2)1^2 /Xv  ~  {SB2  /  B^,)l^2{N/2)1^2 ,  where  N  =  z/A*,  is  the  number  of  wiggler  periods 
the  electron  beam  has  traveled  and  where  zc  ~  A^/2  has  been  assumed.  For  example, 
{SB2 fB„Yt2  =  0.1%  and  N  =  200  gives  {Sx2)if2 /Xv  ~  10-2. 

This  random  transverse  motion  of  the  beam  centroid  leads  to  variations  in  the  parallel 
beam  energy*  through  conservation  of  energy.  The  above  analysis  indicates  that  both  the 


20 


mean  and  the  square  of  (lie  variance  of  (lie  parallel  energy  variation  scale  linearly  with  c, 
with  or  without  transverse  focusing  (i.e.,  transverse  focusing  does  not  significantly  reduce 
the  parallel  energy  deviation).  The  expression  for  the  parallel  beam  energy  for  a  single 
wiggier  realization  is  given  by  Eq.  (13)  along  with  the  appropriate  expressions  for  60  (such 
as  Eq.  (22)  for  the  case  of  a  helical  wiggier  with  transverse  focusing).  Expressions  for  the 
ensemble  averages  of  the  parallel  energy  are  given  by  Eqs.  (26)  and  (27)  for  the  helical 
case.  In  the  limit  (2kpz)2  >>  1,  then  the  normalized  variance  of  the  parallel  beam  energy 
<7||  =  ( (7jj )  —  (7||)2)1/2/7||o  ca"  be  written  as  <7||  ~  ir(6B2 /I?^)1/2 N1^2 ,  where  o2,  >>  1  has 
been  assumed.  For  example,  (SB2 / B^)1^2  =  0.1%  and  N  =  100  gives  d||  ~  3  x  10~2. 

These  effects  may  degrade  FEL  performance.  £or  example,  the  random  walk  of  the 
electron  beam  centroid  may  become  too  large  over  a  sufficiently  short  distance  such  that 
the  radiation  beam  is  no  longer  optically  guided.7  When  this  occurs,  the  electron  beam 
centroid  no  longer  overlaps  the  centroid  of  the  radiation  beam,  thus  leading  to  a  loss  of  the 
FEL  interaction  and  loss  of  FEL  gain.  Ideally,  it  may  be  desirable  to  keep  the  magnitude  of 
the  transverse  beam  displacement  less  than  the  beam  i\ldius  (^xj/rj,  <  1  in  order  to  avoid 
such  detrimental  effects.  Alternatively,  the  parallel  energy  variation  induced  by  the  wiggier 
field  errors  may  become  sufficiently  large  so  as  to  destroy  the  FEL  resonant  interaction. 
Statistically  speaking,  one  can  interpret  the  parallel  energy  variation  due  to  field  errors 
as  an  effective  energy  spread  which  increases  with  increasing  axial  distance.  As  discussed 
above,  c\\  ~  z2/2  and  for  rms  random  field  errors  of  0.1%,  then  after  100  wiggier  periods 
there  exists  an  effective  energy  spread  of  <T||  ~  1.0%,  which  is  a  significant  amount.  In 
order  to  avoid  loss  of  the  FEL  resonance,  it  is  necessary  for  the  parallel  energy  spread  to 
be  small  compared  to  the  intrinsic  power  efficiency  -q  for  FELs  in  the  low  gain  or  high 
gain  regimes,  d||  <  tj.  In  the  trapped  particle  regime,  the  parallel  energy  spread  needs  to 
be  small  compared  to  the  width  of  the  FEL  ponderomotive  potential9  $p,  which  implies 
d||  <  |e$p|/(7m0c2).  For  example,  in  the  low  gain  regime  tj  =  1/(2 N)  which  implies  that 
the  rms  field  errors  must  be  less  than  (SB2 / B\,)1^2  <  l/(2nNs^2),  where  the  expression 
for  <7||  given  in  the  previous  paragraph  has  been  used.  For  N  =  100,  this  implies  that  in 
the  low  gain  regime  (SB2 /B^,)1^2  <  2  x  10“*,  which  may  be  difficult  to  achieve  in  practice. 

In  an  actual  wiggier  magnet,  however,  the  field  errors  SB(z)  are  not  entirely  random 
functions  of  :.  In  practice,  once  the  field  errors  from  each  individual  magnet  pole  have  been 
measured,  one  may  be  clever  in  how  these  poles  are  then  arranged  during  the  assembly 
of  the  wiggier  such  that  the  detrimental  effects  of  these  field  errors  are  minimized.15  For 
example,  if  the  field  error  from  a  given  pole  pair  tends  to  deflect  the  electron  beam  in  a 
given  direction,  then  the  next  pole  pair  should  be  chosen  such  that  it  deflects  the  electron 
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beam  in  the  opposite  direction  so  as  to  keep  the  beam  as  close  to  axis  as  possible.  A  typical 
figure  of  merit  used  in  practice  is  the  line  integral  of  the  field  errors  |  J  dz'6B{z')\  which  is  to 
be  minimized  during  wiggler  assembly.  ]f  the  magnitude  of  the  centroid  displacement  is  to 
be  minimized,  however,  then  the  above  analysis  indicates  that  the  pole  arrangement  should 
be  chosen  to  minimize  the  integral  |  J  dz1  sin  kp{z'  —  z)6J5(r')j,  as  is  shown  by  Eq.  (21).  In 
some  FEL  applications  it  may  be  the  case  that  the  reduction  in  gain  resulting  from  finite 
field  errors  is  dominated  by  the  parallel  energy  deviation  as  opposed  to  the  random  walk  of 
the  electron  beam  centroid.  In  such  cases  it  may  be  desirable  to  arrange  the  magnet  poles 
such  that  the  expression  for  the  magnitude  of  the  parallel  energy  deviation  is  minimized, 
|/>7|||  =  |7||  —  7||ol)  where  7||  is  given  by  Eq.  (13)  along  with  the  appropriate  expressions  for 
bflx-  Preliminary  analysis  by  the  authors  suggests  that  if  one  wishes  to  maximize  the  FEL 
gain,  then  the  magnet  poles  should  be  arranged  such  that  the  magnitude  of  the  deviation 
in  the  relative  phase  of  the  electrons  in  the  ponderomotive  wave  |6i/>|  is  minimized,  where 
Sip  ~  J  dz'S^^.  Notice  that  the  centroid  displacement  scales  as  6x  ~  J  dz'Sfix  ,  whereas 
the  relative  phase  deviation  6ip  depends  on  terms  of  the  form  /  dz'(3x o?>fi±  and  of  the  form 
J  dz'6/3±.  Hence,  minimization  of  Sx  does  not  necessarily  correspond  to  minimization  of 
of  6ip. 

External  steering  coils  may  also  be  used  to  reduce  the  detrimental  effects  of  field 
errors.*-5  For  example,  steering  coils  may  be  used  to  periodically  steer  the  electron  beam 
back  on  axis  and  thus  prevent  the  beam  centroid  displacement  from  becoming  too  large. 
Although  the  random  walk  of  the  centroid  may  be  greatly  reduced  in  this  way,  it  is  not  clear 
that  this  will  also  greatly  reduce  the  parallel  energy  variation  67(1  or,  more  importantly,  the 
relative  phase  deviation  Sxp.  Preliminary  analysis  by  the  authors  indicates  that  if  external 
coils  are  used  to  steer  the  eleciron  beam  back  on  axis  after  a  given  distance  1,  such  that 
6x(l)  =  0,  then  the  mean  value  of  the  relative  phase  deviation  {6rp(l))  is  only  reduced  by 
a  factor  of  three,  i.e.,  1/3  times  the  value  in  the  absence  of  external  beam  steering.  A 
more  complete  analysis  of  the  effects  of  field  errors  on  FEL  gain,  including  the  effects  of 
beam  steering,  is  currently  being  pursued  by  the  authors  and  will  be  the  subject  of  future 
publications. 

Acknowledgements 

The  authors  would  like  to  acknowledge  useful  discussions  with  P.  Sprangle.  This  work 
was  supported  by  the  Office  of  Naval  Research,  Contract  No.  N00014-87-f-0066,  through 
the  National  Institute  of  Standards  and  Technology;  and  by  the  U.S.  Department  of  Energy. 


22 


References 

1)  See,  for  example,  Free  Electron  Lasers,  Proceedings  o}  the  91h  International  FEL  Con¬ 
ference,  ed.  by  P.  Sprangle,  C.M.  Tang  and  J.  Walsh  (North-Holland,  Amsterdam, 
1988). 

2)  K.E.  Robinson,  D.C.  Quimby,  J.M.  Slater,  T.L.  Churchill  and  A.S.  Valla,  in  Free 
Electron  Lasers ,  Proceedings  of  the  8th  International  FEL  Conference ,  ed.  by  M.W. 
Poole  (North-Holland,  Amsterdam,  1987),  p.  62;  K.E.  Robinson,  D.C.  Qubnby  and 
J.M.  Slater,  IEEE  J.  Quantum  Electron.  QE-23,  1497  (1987). 

3)  B.M.  Kincaid,  J.  Opt.  Soc.  Am.  B  2,  1294  (1985). 

4)  C.J.  Elliott  and  B.D.  McVey,  in  World  Scientific  Proceedings  of  Undulator  Magnets 
for  Synchrotron  Radiation  and  Free  Electron  Lasers ,  Trieste,  Italy,  June  1987. 

5)  H.D.  Sliay  and  E.T.  Scharlemann,  in  Free  Electron  Lasers,  Proceedings  of  the  l'<h 
International  FEL  Conference ,  ed.  by  P.  Sprangle,  C.M.  Tang  and  J.  Walsh  (North- 
Holland,  Amsterdam,  1988),  p.  601. 

6)  C.M.  Tang,  E.  Esarey,  W.  Marable  and  P.  Sprangle,  presented  at  the  9th  International 
FEL  Conference,  Williamsburg,  VA,  Sept.  14-16,  1987;  also  presented  at  the  Work¬ 
shop  for  Advanced  Wigglers,  Los  Alamos,  NM,  Dec.  3,  1987;  E.  Esarey,  W.  Marable, 
C.M.  Tang  and  P.  Sprangle,  Bull.  Amer.  Phys.  Soc.  33  (4),  1066  (1988). 

7)  P.  Sprangle,  A.  Ting  and  C.M.  Tang,  Phys.  Rev.  Lett.  59,  202  (1987);  Phys.  Rev. 
A  36,  2773  (1987);  in  Free  Electron  Lasers,  Proceedings  of  the  8th  International  FEL 
Conference,  ed.  by  M.W.  Poole  (North-Holland,  Amsterdam,  1987),  p.  136. 

8)  P.  Sprangle,  R.A.  Smith  and  V.L.  Granatstein,  in  Infrared  and  Millimeter  Waves ,  ed. 
by  K.  Button  (Academic,  New  York,  1979),  Vol.  I,  p.  279;  P.  Sprangle,  C.M.  Tang 
and  C.W.  Roberson,  Nucl.  Instrum.  Methods  A  239,  1  (1985). 

9)  P.  Sprangle  and  C.M.  Tang,  Appl.  Phys.  Lett.  39,  677  (1981);  C.M.  Tang  and  P. 
Sprangle,  in  Free-Electron  Generation  of  Coherent  Radiation,  Physics  of  Quantum 
Electronics,  ed.  by  S.F.  Jacobs,  G.T.  Moore,  H.S.  Pilloff,  M.  Sargent  III,  M.O.  Scully 
and  R.  Spitzer  ( Addison-Wesley,  Reading,  MA,  1982),  Vol.  9,  p.  627. 

10)  W.  Marable,  C.M.  Tang  and  P.  Sprangle,  Bull.  Amer.  Phys.  Soc.  32  (9),  1794  (1987). 


23 


11)  See,  for  example,  Basic  Principles  of  Plasma  Physics:  A  Statistical  Approach ,  S. 
ldiimaru  (Benjamin-Cummings,  Reading,  MA,  1973),  p.  237. 

12)  See,  for  example,  Stochastic  Processes  in  Physics  and  Chemistry ,  N.G.  Van  Kampen 
(North-Holland,  Amsterdam,  1981),  p.  223. 

13)  P.  Diament,  Pliys.  Rev.  A  23,  2737  (1981);  T.I.  Smith  and  J.M.J.  Madey,  Appl. 
Phys.  B  27, 195  (1982). 

14)  W.M.  Fawley,  D.  Prosnitz  and  E.T.  Scharlemann,  Phys.  Rev.  A  30,  2472  (1987). 

15)  M.S.  Curtin,  A.  Bhowmik,  W.A.  McMullin,  S.V.  Benson,  J.M.J.  Madey,  B.A.  Riclnnan 
and  L.  Vintro,  in  Free  Electron  Lasers,  Proceedings  of  the  9th  International  FEL 
Conference,  ed.  by  P.  Sprangle,  C.M.  Tang  and  J.  Walsh  (North-Holland,  Amsterdam, 
1988),  p.  91. 


24 


FJCi.  1.  Schematic  of  random  homogeneous  field  errors  />f?(r)  which  have  a  zero  mean, 
finite  variance  (SB2),//2  and  a  finite  autocorrelation  distance  cc  ~  Xw/2. 
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FIG.  3.  A  typical  transverse  orbit  occurring  in  a  single  wiggler  realization  without  trans¬ 
verse  focusing  for  the  parameters  7  =  270,  flyi  —  2,  Ayi  —  5  cm  and  (6Bi)l'2/Bw  =  0.3%. 
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FIG.  4.  A  typical  transverse  orbit  occurring  in  a  single  wiggler  realization  with  transverse 
focusing  for  the  parameters  7  =  270,  au,  =  2,  Au,  =  5  cm  and  (SB^,)1^2  /  Bw  =  0.3%. 
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FIG.  7.  The  rms  transverse  displacement  of  the  beam  centroid  as  obtained  from  simulation 
(solid  curve)  and  from  theorj'  (dashed  curve)  for  the  case  without  transverse  focusing 
(7  =-  270,  au,  =  2,  \w  =  5  cm  and  {8B2w)lf2 /Bw  =  0.3%). 
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<5x2  >  1/2  (cm) 


FIG.  8.  The  rms  transverse  displacement  of  the  beam  centroid  as  obtained  from  simulation 
(solid  curve)  and  from  theory  (dashed  curve)  for  the  case  with  transverse  focusing  (7  = 
270,  aw  =  2,  \w  =  5  cm  and  (SB^)1  / Bu.  =  0.3%). 
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FIG.  10.  The  normalized  mean  parallel  energy  ("y || ) /*7 jto  —  1  as  obtained  from  simulation 
(solid  curve)  and  from  theory  (dashed  curve)  for  the  case  with  transverse  focusing  (7  = 
270,  a-u,  =  2,  Xw  =  5  cm  and  (SBl,)1^2 / Bu,  =  0.3%). 
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FIG.  11.  The  normalized  variance  of  the  parallel  energy  ((-,jj)  -  ("yn)'  )1/2/*)||o  as  obtained 
from  simulation  (solid  curve)  and  from  theory  (dashed  curve)  for  the  case  without  trans¬ 
verse  focusing  (7  =  270,  au.  =  2,  Au.  =  5  cm  and  (tBl,)1'2  Bu.  =  0.37r). 
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FIG.  12.  The  normalized  variance  of  the  parallel  energy  {{‘l;  -  'Oil)2  )1/2/7|io  as  obtained 
from  simulation  (solid  curve)  and  from  theory  (dashed  curve)  for  the  case  with  transverse 
focusing  (->  =  270,  a„.  =  2,  Xu.  =  5  cm  and  {f< D^,)1  / B w  =  0.3%). 
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